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Abstract 



A static spherically symmetric asymptotically flat spacetime may allow for cir- 
cular closed null-geodesics which are said to belong to a photon sphere. In the 
context of gravitational lensing in the strong deflection regime, the presence of a 
photon sphere leads to unbounded angle of deflection of light (multiple turns) and 
bi)' formation of relativistic images. In this paper we show that photon spheres may 

form in some configurations of boson stars constructed with free massive complex 
scalar field nonminimally coupled to gravity. Assuming the boson star is transpar- 
' ent to light, photon spheres would not only give raise to phenomena in the realm of 

\^ • strong gravitational lensing, but also to considerably increased photon flux in the 

! central region of the star, relative to the flux in its surroundings. 

^ : 1 Introduction 



Deflection of light by a massive body, as predicted by General Relativity, led to one of the 
flrst experimental veriflcations of the theory [1]. Deflection of light is also the mechanism 
behind the rich phenomenology of gravitational lensing that has seen a wide range of 
applications in astrophysics [2j. In observationally relevant circumstances gravitational 
lensing occurs in the so called weak deflection regime; the light passes by a massive body 
at distances that are much greater than the Schwarzschild radius corresponding to the 
mass, and the overall angle of deflection is small. However, despite of known difficulties 
on the observational side, there is also considerable interest in gravitational lensing in 
the strong deflection regime, see e.g. [3] and references therein. This primarily includes 
the case of lensing by the Schwarzschild black hole [U El E] , by the regular black hole [7] , 
or other types of black holes, see [S] for a review. Light passing sufficiently close to the 
Schwarzschild black hole may, due to the presence of the photon sphere in the spacetime, 
suffer arbitrarily large deflection as it makes one or more complete turns around it before 
escaping to spatial infinity. In the context of gravitational lensing this leads to the 
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formation of the so called relativistic images. Mathematically rigorous treatment of the 
concept of the photon sphere can be found in [9] . 

In this paper we look for a simple matter model that can be used to construct static, 
spherically symmetric self-gravitating bodies, that curve the spacetime in such a way that 
photon spheres are formed. In this regard one is faced with two possible routes. One 
is to construct a compact body, i.e. a body that is finite in the radial extent, with the 
surface lying below the photon sphere of the Schwarzschild spacetime. If such body is 
assumed to be non-transparent to photons, then its properties as a gravitational lens are 
equivalent to those of a Schwarzschild black hole. The other, more interesting route, is 
to allow the energy-momentum content to take up all space and assume that the matter 
it describes is transparent to light. The simplest matter model for which it is natural 
to assume that it does not interact electromagnetically is the scalar field. For example, 
lensing by configurations of massless scalar field with a naked spacetime singularity was 
studied in [lOj. As we intend to show, photon spheres can form within the regular objects 
formed of scalar fields known as boson stars. 

Boson stars were originally conceived as Klein-Gordon geons [TT1[T2], the self-gravitating 
configurations of the free massive scalar field minimally coupled to gravity, while among 
the extensions of the original model one finds quartic self- interaction [13], non-minimal 
coupling to gravity [HI [15], addition of the Brans-Dicke field [16], etc. Due to their sim- 
plicity boson stars proved to be an attractive theoretical model for studying properties 
of self gravitating objects in General Relativity. They were also considered as candidates 
for astrophysical objects such as supermassive Sgr A* at the centre of our galaxy [T7] . 
Several reports on the subject appeared over the past decades [IHI [El EO] . Gravitational 
lensing by boson stars was considered in [211 122]; relatively large angles of deflection of 
light forming primary and secondary images were found, but the possibility of formation 
of photon spheres and relativistic images was not reported so far. 

The paper is organized as follows: In Section [2] we give a brief overview of the prop- 
erties of null-geodesics in a static, spherically symmetric, asymptotically flat spacetime 
with photon spheres, and we introduce the notation that will be used through the rest of 
the paper. We also discuss a simple analytical model of a compact body with the surface 
within the Schwarzschild photon sphere. In Section [3] we construct boson stars from free 
massive complex scalar field nonminimally coupled to gravity and we show that photon 
spheres form some configurations. We sum up in Section [H 

2 Notation and a simple model 

We will be considering the behaviour of null-geodesics in static, spherically symmetric, 
asymptotically fiat spacetimes, free of singularities or black holes. The line element can 
be written using the coordinates x"" = [t, r, i), ip) and geometrized units (c = 1 = G) as 

ds^ = gab dx'^ dx' = -e2*(^)dt2 + e^^^^Mr^ + rW, (1) 

where the radial coordinate r is the area radius and dfi^ = dt?^ + sin^ dip"^ is the line 
element on the unit sphere. The Schwarzschild spacetime contains a black hole of mass 
M > with the event horizon at r = 2M. Outside the black hole the metric components 
of the Schwarzschild spacetime are given by 

e2*(r) ^ g-2A(r) ^ ^ _ ^ ^ 2M, (2) 
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M being the ADM mass of the spacetime. We will encounter the situation where the 
the energy-momentum tensor vanishes at radii greater than some constant R > 2M and 
where the metric components are given by (I2j) exactly, as well as situations where nonzero 
energy-momentum takes up all space and where ([2]) only represents the asymptotic form 
of the metric as r — )■ cxd. 

The Lagrangian for the geodesies constrained to the -d = n /2 plane can be written as 



where the overdot denotes differentiation with respect to the parameter of the geodesic 
and the choice of the sign depends on whether the geodesic is spacelike (+) or timelike 
(— ). As L does not depend on t and there are two constants of motion, dL/dt and 
dLjdup. The ratio of the two constants, 

6 = ^ = -rV-«^ = const., (4) 
dLjdt dt ' ^ ^ 

is proportional to the angular momentum of the test particle following a timelike geodesic 
and is known as the impact parameter. For a geodesic that can be extended to spatial 
infinity, i.e. to the region where the spacetime is essentially flat, the absolute value of 
the impact parameter can be interpreted as the distance between the considered geodesic 
and the purely radial geodesic that is (at spatial infinity) parallel to it. 

By continuity, the impact parameter h defined in is a constant of motion also for 
the null-geodesics, where in addition we have L = 0. Eliminating d(/9/dt one obtains an 
equation for dr/dt that can be written in the form of the nonrelativistic energy equation, 

1 /dr^ ^ 
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where £" = and the 'effective potential' V(\)^ r) is given by 

^(^-) = ^^ l)- (6) 

In analogy with the familiar scenario from classical mechanics we see that null-geodesics 
exist only where 1^(6, r) < 0, as the r.h.s. of ([5]) must be non-negative for dr/dt to be 
real. If V^(6, tq) = then at r = tq there is a turning point of r(t) if dV{\)^ r)/dr\r=ro 7^ 0, 
or an equilibrium point if dV{h, r)/^r\r=J.^^ = 0. A turning point is an inner turning point 
(a minimum of r(t)) if dV(h,r)/dr\r=TQ < 0; it is an outer turning point (a maximum 
of r{t)) if dV{b,r)/dr\r=ro > 0. An equilibrium point of r(t) implies existence of closed 
(circular) null-geodesics along equators of the hypersurface r = tq which is known as the 
photon sphere. A photon sphere is said to be stable if d'^V{b,r)/dr'^ > 0, or unstable if 
d'^V{b,r)/dr'^ < 0. In the neighborhood of a stable photon sphere there are bound null 
geodesies with oscillatory r(t) such that the outer turning point is above, and the inner 
turning point is below the photon sphere. Apart from circular null-geodesics, an unstable 
photon sphere implies existence of null-geodesics that asymptotically approach the photon 
sphere, either from the inner or from the outer side, making an infinite number of turns 
around it. It can be shown that the outermost photon sphere in an asymptotically flat 
static spherically symmetric spacetime is an unstable photon sphere. 
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To access the properties of nuU-geodesics in a given static spherically symmetric space- 
time in a systematic way it is convenient to consider the relation among the impact pa- 
rameter b and the radial coordinates of turning points or photon spheres tq. This relation 
follows from the condition dr/dt\r=ro = 0, or V{b,ro) = 0, and reads 

b{ro) = ±roe-*(^°). (7) 

The double sign in ([7]) corresponds to the sign of the angular momentum, or the time 
direction one chooses along the given null-geodesic. For simplicity of notation, with no 
loss in generality, we will assume positive b in what follows. It can be easily shown 
that solutions to ([7j) with d6/dro > represent inner turning points, while those with 
d6/dro < represent outer turning points of null-geodesics; the solutions with dfe/dro = 
and d^6/drQ < represent stable, while those d^6/drQ > represent unstable photon 
spheres in the spacetime. 

Let us first consider the region of the Schwarzschild spacetime outside the black hole 
where the metric components are given by ([2]). There is no solution to ([7j) with ro > 
for impact parameter b < S^^^M. This implies that null-geodesics with b < S'^^^M extend 
to spatial infinity at one end, while their other end reaches the black hole horizon. At 
b = S'^/^M there is the unstable photon sphere with area radius tq = 3M (see point A 
in Fig. [T]). For b > one finds two turning points. The inner turning point with 

ro > 3M and with the asymptote ro = b — Masb-^oo represents the closest approach 
of the null-geodesics extending to spatial infinity. The outer turning point with ro < 3M 
and with the asymptote ro = 2M as 6 — )■ oo represents the outermost point reached by 
the geodesies with both ends at the black hole horizon. 

We now proceed to replace the central segment of the Schwarzschild spacetime con- 
taining the black hole of mass M with an exact solution of the Einstein equations repre- 
senting a static spherical body (compact object) of mass M and surface radius R > 2M. 
We choose to do so in such a way that the joining hypersurface r = R lies inside the 
Schwarzschild photon sphere, i.e. we require 2M < R < 3M. This implies that the ratio 
2M/R, known as the surface compactness of the body, belongs to the highly relativistic 
regime, 

2 2M 

Here one recalls Buchdahl's upper bound on the surface compactness for bodies involving 
isotropic pressures in their interior [23], 2M/R < 8/9, and also that many stellar models 
develop disturbing features such as violation of the energy conditions or loss of dynamical 
stability at high values of surface compactness. However, simple models of massive bodies 
with sufficiently high surface compactness to satisfy the condition ([8]) can be found and 
we choose here to work with a model consisting of a massive spherical shell of infinitesimal 
thickness supported against gravity by its own surface pressure. It was shown in Ref. [21] 
that for the case of the flat metric in the shell interior this model allows the surface 
compactness as high as 24/25 without the violation of the dominant energy condition. 
The metric components inside the shell are given by 



e 



^^CO = 1 _ 2M/R, e'^^^) = 1, r<R. (9) 



It is important to note that joining the metrics ([2]) and ([9]) at r = i? is not smooth. While 
gtt = — e^**^''^ is continuous, grr = e'^^^^^ has a discontinuity at r = R. This feature is 
closely related to the 5-shaped distribution of matter of the massive shell along the radial 
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coordinate. The full technical details of the construction of the energy-momentum tensor 
on the hypersurface r = R which involves the application of Israel's thin shell formalism 
[25] are not needed for the present discussion; we refer the interested reader to |23j or to 
f2E\ [27] where a similar model with a segment of the de Sitter spacetime in the interior, 
known as the gravastar, is worked out in detail. 

The relation among the impact parameter b and the radial coordinate of the turning 
point or a photon sphere tq for the spacetime containing the massive shell obeying (j8]) 
is shown in Fig. [1] The curved line is the plot of ([7]) for the Schwarzschild metric (|2]) 
which is relevant for rg > -R (the part of this curve below the point B is shown dashed). 
The straight line with one end at 6 = tq = and the other at B is the plot of ([7]) 
for the metric inside the shell which is relevant for tq < R. The point A with the 
coordinates tq = 3M and b = is the Schwarzschild photon sphere, while the part 

of the curve above A represents the inner turning points (radii of closest approach) for 
null-geodesics with b > reaching spatial infinity. These nuU-geodesics are entirely 

contained in the Schwarzschild segment of the spacetime and are in no way affected by 
the massive shell replacing the black hole. Part of the curve between A and B represents 
the outer turning points, and the part of the straight line between C and B represents 
the inner turning points, for bound null geodesies with impact parameter b in the range 
< b < bs = R/ \/^1^-2^MJR. Part of the straight line below C represents the inner 
turning points (radii of closest approach) for null-geodesics reaching spatial infinity with 
b < 3^/2M. The turning points of these null-geodesics lie within the sphere of area radius 

rc = - 2M/R, (10) 

which is smaller than the radius of the shell. It is easy to see that rc — )■ for 2M/R — )■ 1. 
This also means that there is no null-geodesic extending to spatial infinity with the radius 
of closest approach tq in the range rc < < 3M. The behaviour of a bundle of null- 
geodesics that are parallel at spatial infinity is shown in Fig. [2] as it reaches the massive 
shell of radius R = 5M/2. The null-geodesics with the impact parameter b < (solid 
lines) show strong degree of focusing in the interior of the massive shell. This interesting 
feature could not be confirmed analytically since the null-geodesics in the region r > R 
could only be constructed by numeric integration of (jlj) and ([5j). 

In summary of the model we have considered, we refer to the points A, B and C 
in Fig. [TJ Point A is the outer (unstable) photon sphere, while B is the inner (stable) 
photon sphere. Point C is defined so that be = bA- The null geodesies extending to spatial 
infinity have turning points with rg < rc for impact parameter b < b^, or with Tq > 
for b > b^- In addition, there is a class of bound null geodesies with impact parameter b 
in the range bA < b < bs and the turning points tq in the range rc < r^ < rA- All these 
qualitative features, i.e. the points A, B and C arranged as in Fig. [H are to be found in 
a different context in the following Section. 



3 Photon spheres in boson stars 

The action involving the massive complex scalar field (p nonminimally coupled to gravity 
is 

S = J d'xy^ + ^0*0^ (V>*)(V„0) -/xV^) , (11) 

where ^ is the nonminimal coupling constant, R is the curvature scalar and /i is the field 
mass. Since the action is invariant under the global phase transformation — )■ e"^0, there 



5 




M 2M 3^I^M 
impact parameter b 



Figure 1: Turning point or photon sphere radius tq vs. impact parameter h (thick sohd 
hne) for null-geodesics in the spacetime containing the thin spherical shell of mass M 
and radius R = |M. (Null-geodesics exist only to the left of the thick solid line; shaded 
region corresponds to the shell interior.) Thick dashed line is the continuation of the 
tq vs. b curve for the case of black hole of mass M. Thin dashed lines are provided as 
eye-guides, the line tq = b — M is the asymptote of 6(ro) for 6 — )■ oo. 




Figure 2: A bundle of asymptotically parallel null-geodesics passing by (impact parameter 
b > dashed lines) or through (6 < solid lines) the spherical shell of mass 

M and radius R = |M (thin solid circle with shaded interior). The turning point of each 
null-geodesic is indicated with the small circle. 
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is the conserved current, ja = i((Va</>*)0 — (Va0)0*), and the corresponding generator, 
the particle number A^. The variation of f lTT]) with respect to the metric gives the Einstein 
equation. Gab = ^T^Tab, where Gab = Rab — \gabR is the Einstein tensor and the energy- 
momentum tensor is given by 

Tab = V,0* Vb0 + Va0Vb0* - gab (V^^^V^^ + 

- 2^ (l)*(l)Gab - 2i gabWc(l)*(l) + 2^ VaVb0*0. (12) 

Note that the energy-momentum tensor can be written in a different way if the term 
2C,(f)*(f)Gab is transferred to the l.h.s. of the Einstein equation, which is then divided by 
l + 167r^0*0. There are also different approaches in the in the interpretation of the energy- 
momentum tensor in case of the nonminimally coupled scalar field, for a discussion see 

[2H]- 

In spherical symmetry, with the metric written as in ([1]), the nontrivial components 
of the Einstein tensor are 

G\ = (e-2^(l - 2rA') - l) , (13) 
G\ = r-' (e-2^(l + 2r$') - l) , (14) 
6-% = G^^ = r-2e-2^((r$' - rA')(l + r$') + r^^") (15) 

(explicite notation of the r-dependence of functions is omitted and the prime denotes 
the r-derivatives). With the stationary Ansatz for the complex field, = 0(r)e~"^*, 
where from here on denotes the real function, the nontrivial components of the energy- 
momentum tensor can be identified with the energy density p, radial pressure p and the 
transverse pressure q. They are given by 

p = -T\ = (/i2 + e-2*u;2)^2 ^ g-2A^,2 ^ 2^^2^*, 
2A / 1/2 



+ 4ee-^^(0'^ + (2/r - A')0> + <p"<p), (16) 
p = T\ = {-p^ + e-2*a;2)^2 ^ g-2A^/2 _ 2^02^.^ 

-4ee-2^(2/r + <l>')0V, (17) 
q = T\ = T\ = (-/i^ + e-2*a;2)^2 _ ^-2A^n _ 2ict>^G% 

- 4ee-2^ (0'' + (1/r - A' + $')0'0 + 0>) . (18) 

Plugging f|T3l) - f|T5l) and f|T6|) - f|T8|) into the Einstein equation gives the system of three 
coupled ordinary differential equations with the three unknown functions: $, A and 0. 
The differential order of the system equals five since the highest derivatives that occur 
are A' and 0". The conservation condition VaT"?, = 0, which gives 

r'^e^^i-^? + e-^*a;^)02 + (2 - rA' + r<l>')r0'0 + rV'0 

= 2^(1 - e^^ + 2r$' + r^^'^ - rA'(2 + r<l>') + r2$")02, (19) 

can be used to eliminate $" from the system. This reduces the differential order of the 
system by one. (Note that (|T9l) is equivalent to what one obtains from (fTTj) for the 
equation of motion for the scalar field.) The solution is determined with the boundary 
conditions which reflect the requirements that all functions are regular at r = 0, that 
the metric is asymptotically fiat and coincides with ([2]), and that the field vanishes as 
r — )■ oo. The number of the required boundary conditions equals the differential order 
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of the system (four) plus one since the system also involves the unknown frequency uj 
(eigenvalue) which has the role of an additional degree of freedom. The particle number 
is given with 



and can be obtained once the functions $, A and 0, and the frequency oj have been 
determined. 

In order to set up the numerical procedure several further steps are taken. It is 
convenient to replace the metric function A(r) with the so called 'mass function' m(r) 
defined with 



(in the non-relativistic regime m is the mass within the sphere of radius r). The asymp- 
totic value of m as r — )■ oo is the ADM mass of the spacetime, M. Radial coordinate 
X = ixr is introduced which is then mapped onto the compact domain with the transfor- 
mation X = x/{l + x) G [0, 1]. Following the usual practice in the literature we introduce 
the field variable a = y/Sncj). Finally, for the boundary conditions at x = (correspond- 
ing to r = 0) we choose m = 0, a = ao and a' = 0, while at a; = 1 (corresponding to 
r — oo) we choose $ = and a = 0. The central value of the field variable, ctq = (j{0), 
is used to parametrize the family of solutions obtained for some ^; for each ^ and ctq, the 
solution consists of the three functions, $, m and a, and the eigenvalue u. We used the 
COLSYS boundary value code [29] and checked that our procedure fully reproduces the 
results presented in the Table 1 of Ref. pi ]. 

We limited our search for the photon spheres in the spacetimes sourced by the boson 
stars to the configurations without nodes in 0. (The configurations with one or more 
nodes are sometimes referred to as the excited boson stars [SO]-) For a given value of ^, 
as the central value of the field variable ao increases starting from ao = 0, the mass M 
and the particle number N also increase from zero up to the (first) critical configuration 
where they reach their maxima. In analogy with the similar behaviour of mass in fluid 
stars, the first maximum of M suggests the onset of dynamical instability. In the case 
of minimal coupling of the scalar field and gravity = 0) the loss of stability with 
respect to radial perturbations at the first critical configuration has been confirmed by 
perturbative analysis [311 [32] cilso by other methods [33l[3l]. As do increases beyond 
the first critical configuration, the mass M and the particle number go through a 
sequence of minima and maxima resembling strongly damped oscillations about respective 
asymptotic values. Subsequent maxima of M and A^, that will also be referred to as 
critical configurations, are lower than their maxima at the first critical configuration. 

The binding energy of a boson star can be defined as 



For all values of ^ that we considered we found that the critical configurations (maxima 
of M and A^) correspond to the minima of E'b. In all first critical configurations the 
binding energy E\y is negative which means that the bosons in such configurations may 
not disperse to infinity without the input of additional energy. For negative ^ and up 
to ^ = 2, as (To increases beyond the first critical configuration, Ey, becomes positive 
before the second critical configuration is reached and remains positive with the further 
increase of (Tq. Interestingly, for ^ = 4 and greater, at the second critical configuration 
El, is negative. The dependence of the binding energy (122|) on ctq for some values of C, 




(20) 



Eb = M - fiN. 



(22) 
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is shown in Fig. |3l one can see that for ^ = 4 there is a positive binding energy barrier 
between the first and the second critical configuration, while for sufficiently large values 
of ^ this barrier disappears. It should be emphasized, however, that negative E-^ at some 
equilibrium configuration is a common characteristic of a gravitationally bound states, 
but does not by itself imply stability. 

As outlined in the preceding section, in a given spacetime the photon spheres can 
be located by examining the behaviour of the function 6(ro) defined in ([7]). If &(ro) 
is monotonically increasing with ro > 0, there is no photon sphere in the spacetime, 
whereas if extrema exist, they correspond to photon spheres. In the present context a 
more robust procedure for locating the photon spheres can be formulated. The condition 
for the extremum of 6(ro) can be written as 

= ^ = e-*(l-r$'), (23) 
dro 

where $' can be expressed in terms of the mass function m and the radial pressure p 
using the well known equation 

= i — 

r2(l — 2m/r) 

(in the non-relativistic regime $ is the gravitational potential). It follows that the con- 
dition for the existence of the photon sphere in the spacetime can be expressed as 

— = -il-Anr^p). 25 
r 6 

The quantity on the l.h.s. is known as the 'compactness function' and is always less 
than unity, while the r.h.s. is less than two-thirds if the radial pressure is positive. It 
is interesting to observe that when one considers the compact object of mass M and 
surface radius R, the condition (|8]) for the photon sphere in the vacuum segment of the 
spacetime requires the surface compactness 2M/R to be greater than 2/3; here we find 
that, for positive radial pressure p, the compactness function may satisfy the condition 
at values that are less than 2/3. 

In the range of values of the non-minimal coupling constant ^ that we considered 
(see Table H]), we did not find photon spheres at or below the value of ao corresponding 
to the first critical configuration. Interestingly, the photon spheres were found in the 
second critical configurations for ^ = and above. However, only at ^ = 4 and above 
we have negative binding energy at the second critical configuration. Fig. H] shows the 
relation among the turning point radius tq and the impact parameter b for the boson 
star with ^ = 4. Following the notation scheme introduced in the preceding section, the 
outer (unstable) photon sphere is marked with A and the inner (stable) photon sphere is 
marked with B. The impact parameter corresponding to the photon sphere A is denoted 
with bps. The area radius of the sphere containing the turning points of all null-geodesics 
with b < bps is denoted with r^. The layout of the points A, B and C for higher values 
of ^ is very similar to the layout shown in Fig. HJ Fig. [5] shows how the condition (125|) 
can be used to locate the photon spheres. 

In the simple model of the spacetime with photon spheres constructed in the preceding 
section we found strong degree of focusing of null-geodesics in the interior of the compact 
object. Here, as a measure of shrinking of the bundle of null-geodesics that are initially 
parallel at spatial infinity, one can consider the ratio of the two proper areas defined as 
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(24) 




Figure 3: Binding energy ( 122|) of the boson star per particle in units of the boson mass 
/i vs. the central value of the field variable (Tq for the nonminimal coupling constant 
^ = 2,4,8,16. First and second critical configurations are indicated with filled and 
hollow circles, respectively. 



follows. The first is the proper area of the circular cross section of the part of the bundle 
at spatial infinity containing all null geodesies that (in the course of time) pass through 
the photon sphere and enter its interior. These geodesies must have the impact parameter 
b < bps so the proper area of the cross section is A^o = b'^^^. The second proper area is 
that of the sphere containing the turning points of the geodesies that enter the photon 
sphere. As the area-radius of this sphere is r^, its proper area is Aq = AtqIt. Defined in 
this way, the null geodesies that pass through A^o also pass through Ac, and vice versa. 
The ratio of the two proper areas is 



A 



CO 



C ^'c' 



Ar 4:r^ 



(26) 



It can be interpreted as the ratio of the number of null-geodesics passing through unit 
proper area at r = r^; and at r — )■ oo. As it can be seen from the figures in Table [H in 
the configurations of boson stars where we found photon spheres this ratio is by several 
orders of magnitude greater than unity. If one also considered the number of null-particles 
(photons) passing through unit proper area in unit proper time, assuming steady flux at 
r — )■ oo, an additional factor of e"**^'''^-' = bp^/rc would enter due to the difference in clock 
ticking rates. 



4 Conclusions 

We have shown that photon spheres are present in some conflgurations of the boson 
stars constructed with the free massive scalar field nonminimally coupled to gravity. 
This implies the existence of null-geodesics that make arbitrarily many turns around the 
central region of the star before escaping to spatial infinity. Assuming the boson star is 
transparent to light, this further leads to the possibility of formation of the relativistic 
images in the strong deflection regime of gravitational lensing. With the finding of the 
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Table 1: Central value of the field variable a^, mass M and the particle number of 
the first and the second critical configuration of the boson star obtained with a sequence 
of values of the nonminimal coupling constant ^. The impact parameter corresponding 
to the photon sphere bps, area radius of the photon sphere r^, and the area radius rc 
of the sphere containing the turning points for nuU-geodesics with b < bps are given for 
configurations where the photon sphere exists. Mp denotes the Planck mass. 





(Co)crit. 


M/(M|//i) 


iV/(M|//x2) 


bps/M 


ta/M 


rc/M 


16 


0.146287 


2.89786 


3.48375 










0.273464 


2.18529 


2.27896 


3.50852 


0.196796 


0.0233624 


8 


0.203279 


2.05492 


2.44047 










0.383436 


1.54653 


1.59122 


3.47565 


0.197621 


0.0236546 


4 


0.277568 


1.46454 


1.70347 










0.533207 


1.09811 


1.10413 


3.41338 


0.199297 


0.0242506 


2 


0.365541 


1.05839 


1.19175 










0.729797 


0.78796 


0.764328 


3.29866 


0.202721 


0.0254954 





0.271059 


0.633001 


0.653003 










1.62527 


0.382621 


0.330516 


1.74432 


0.176402 


0.134695 




Figure 4: Radius of the turning point (or photon sphere) ro vs. impact parameter b for the 
nuU-geodesics in the spacetime containing the boson star with the nonminimal coupling 
constant ^ = 4 at the first (no photon spheres) and at the second (photon spheres A and 
B) critical configurations. The impact parameter corresponding to the photon sphere 
A is indicated with bps, the area-radius of the sphere containing the turning points for 
nuU-geodesics with b < bps is indicated with rc- 
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Figure 5: Locating of the photon spheres A and B in the spacetime containing the boson 
star with the nonminimal couphng constant ,^ = 4 at the second critical configuration as 
the solutions of fl2B|) . 

photon spheres within the boson stars, the scalar field appears to be a suitable matter 
model for modeling of strong gravitational lenses that do not involve singularities of black 
holes in the spacetime. Another interesting feature of global (no black holes) spherically 
symmetric spacetimes with photon spheres, relative to those with black holes, is the high 
degree of focusing of the incident light into the small central region of the spacetime. For 
the boson stars we have considered we were able to give a measure of contraction of the 
incident bundle of parallel nuU-geodesics in form of the ratio (126|) . 

However, one must not overlook the not so desirable properties of the configurations of 
boson stars in which photon spheres were found. In the case of the minimally coupled field 

= 0) photon spheres were found in configurations that are known to be dynamically 
unstable and that have positive binding energy (positive work is required to 'compress' 
the bosons with initially zero kinetic energy at spatial infinity into the equilibrium config- 
uration). As we have shown, the unusual sign of the binding energy of the configurations 
with photon spheres can be reversed with the introduction of sufficiently strong nonmin- 
imal coupling of the scalar field to gravity, but the issue of instability remains open. In 
this regard it would be interesting to carry out the linearized perturbation analysis for 
the boson stars with nonminimally coupled field (e.g. like it has been done for ^ = in 
[32]) and if unstable modes are found in the configurations with photon spheres, to ob- 
tain the timescales of the growing perturbations. If these are found to be long relative to 
some time scale refiecting the size of the spacetime region in which gravitational lensing 
is taking place, then the boson stars with photon spheres would not be ruled out on the 
basis of instability as models of massive dark objects giving raise to phenomena of strong 
gravitational lensing. The possible extension of this work could be to explicitly study 
the angular pattern, magnifications and the time delays of relativistic images formed by 
boson stars and compare them with the results that are obtained when the massive dark 
objects in the centres of nearby galaxies are modeled as black holes [HI ES]- 
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